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ABSTRACT 



An attempt is made to investigate the relationship be- 
tween the digitized satellite brightness data and the 
kinematically computed 200-rab divergence in the tropical 
western North Pacific during the period 19 April-20 December 
1971, Both spectral and correlation analyses techniques are 
used. The power spectra of the two fields reveal two domi- 
nant period bands: one centered at - 10 days and the other 
at -- 5 days. Horizontal structures, as shown in the inter- 
longitudinal cross-spectra in the two bands, suggest that 
brightness is more organized than divergence, although the 
two fields resemble well each other whenever an organized 
divergence pattern is observed. Inter-parameter cross- 
spectra and correlation results indicate a generally in- 
phase relationship. However, the coherence squares and 
correlation coefficients are lower than those found by 
Wallace (1971) for the 4- to 5-day waves using direct radio- 
sonde data. 
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INTRODUCTION 



I- 

Satellite photographs have revealed that much of the 
cloudiness in the tropics is in the form of large, con- 
nected masses of very bright clouds. These cloud masses 
consist of varied individual cumulus clouds. The extensive 
cirrus canopies observed on top of these cumulus clouds 
account for most of their brightness, and may last a few 
hours longer than the cumulus cells themselves. These 
’'cloud clusters,” as they are now called, have attracted 
a great deal of attention for the past few years since, 
hypothetically, they play an important role in the general 
circulation of the atmosphere and in the formation of 
tropical storms. Chang (1970) and Wallace (1971) observed 
that many of the cloud clusters in the tropical Pacific 
are well organized and have characteristics which resemble 
those of synoptic-scale, westward propagating wave distur- 
bances. Sikdar and Soumi (1971) and Reed and Recker (1971) 
found evidences that these cloud clusters are envelopes 
of deep convection in regions of organized upward vertical 
motion associated with synoptic-scale waves, Sikdar and 
Soumi (1971) further suggested that the deep convection 
cells tend to concentrate in a small area and are generally 
located in the geometrical center of the cloud cluster; 
moreover, they postulated that the time variations of these 



9 



cells could be excellent indicators for the convective 
energy transport to the tropical upper troposphere. 

Williams and Gray (1973) and Yanai et al (1973) also indi- 
cated that the active cloud clusters exhibit vertical 
characteristics similar to those found by Reed and Recker 
(1971) . 

The paucity of meteorological stations in the vast 
tropical oceans accentuates the imperative for further ex- 
ploiting the utilization potential of available satellite 
data as aids in tropical analysis. In particular, because 
of the relationship between cloud cluster and synoptic- 
scale tropical motions, the digitized cloud brightness data^ 
present a prime research area for such a purpose. 

The physics of the atmosphere suggests a close relation- 
ship between digitized satellite cloud brightness (hereafter 
referred to, for brevity, as brightness) data and the large- 
scale vertical motion field. Williams and Gray (1973) 
pointed out that the cumulonimbi do serve as agents in the 
production of the large, bright cirrus canopies of the cloud 
clusters, with the developing and conservative clusters 
generally capable of maintaining their cumulonimbi in the 
presence of a steady low-level mass convergence, otherwise 
they gradually die out. Since low-level moisture over the 
oceanic tropics is readily available, the large-scale upward 



Digitized brightness data have been made available 
since 1967, through the Advanced Vidicon Camera Systems 
of operational satellites. 
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motion provides a favorable environment for the development 
and maintenance of these clusters. Furthermore, the brighter 
clouds are usually thicker and, therefore, may indicate a 
larger amount of latent heat release (Gruber, 1974). Since 
deep tropical systems are characterized by temperature fluc- 
tuations which, by and large, are one order of magnitude 
less than the rate of latent heating, thermodynamic energy 
balance requires that the large-scale vertical motion be 
approximately proportional to heating [Wallace (1971) and 
Holton (1972)]. A jprj.orj., therefore, a positive correlation 
between brightness and large-scale vertical motion may be 
expected by reason of heating alone. 

In fact, indications of good correspondence between area- 
averaged brightness and synoptic-scale vertical motion 

associated with the 4- to 5-day tropical waves in certain 
2 

sections of the tropical western North Pacific has been 
found by Wallace (1971), on which basis he postulated that 
the entire vertical motion field over large tropical areas 
may be estimated from satellite data alone. If this is the 
case, brightness is certainly one of the best available in- 
dicators of tropical weather disturbances considering that, 
at present, there is no direct way of obtaining a reliable 
vertical motion field in the data-sparse tropical regions. 



These are in the Kwa j alein-Eniwe tok-Ponape (KEP) 
triangle of the Marshall Islands and the Guam-Tr uk-Yap (GTY) 
triangle farther west. 
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The correspondence between brightness and vertical mo- 
tion may also be extended to large-scale horizontal diver- 
gence at levels where divergence is large, since several 
observational studies [Wallace (1971), Reed and Recker (1971), 
Nitta (1972), Williams and Gray (1973), and Yanai et al 
(1973)] indicated similar profiles of the vertical motion 
or divergence field over convectively active tropical areas. 
Implicitly, the same conclusion can be drawn from the verti- 
cal profiles shown by Hayashi (1974) from analysis of tropi- 
cal output from the general circulation model of the NOAA 
Geophysical Fluid Dynamics Laboratory at Princeton. This is 
also consistent with the results of the composite study of 
Williams and Gray (1973), which show a clear in-phase rela- 
tionship between cloud clusters and upper-level divergence. 

If satellite data could be used as a crude representation 
of the upper-level divergence, such data may be very useful 
in large-scale analysis over data-sparse tropical regions, 
particularly so when we consider that the upper-level circu- 
lation pattern is largely controlled by the strong divergence 
near the 200-mb level. By use of the quasi-barotropic 
property of the large-scale tropical flow, it may be possible 
to utilize the divergence information that could be inferred 
from the brightness data as a diagnostic tool in examining 
the flow pattern near the strong divergent level in a manner 

similar to Holton and Colton’s (1972) model for seasonal mean 
3 

motions. With infrared (IR) radiation data, supplied by the 
3 

A diagnostic model of using brightness as an aid for 
large-scale tropical flow is discussed in Jacobs (1974). 
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recently developed satellites, the possibility of represent- 
ing the divergence field by brightness has increased. This 
Is seen from the fact that one can determine cloud-top tem- 
peratures from IR data. These cloud-top temperatures, which 
are indicative of the vertical extent of cumulus development, 
can then be used to calculate the resulting latent heat re- 
lease [Koffler et al (1973) and Rao (1970)]. Thus, it is 
possible to modify the brightness field and improve the esti- 
mate of the upper-level divergence (and vertical motion) 
field. Motivated by these considerations, this study seeks 
to extend the investigation of Wallace (1971) on the relation- 
ship between brightness data and upper-level divergence over 
a broader geographical area. In this study, we hope to find 
out the extent to which the current brightness data can be 
used to represent the synoptic-scale divergence in the tropi- 
cal upper troposphere. In so doing, both spectral and corre- 
lation analyses are used. 

The main reason for performing spectral analysis is be- 
cause of the importance of synoptic-scale tropical distur- 
bances which may highly influence the variations in both the 
divergence and brightness fields. These variations may be 
characterized by certain dominant time scales associated 
with synoptic-scale disturbances which may or may not be 
related to the easterly waves. It is the interest of this 
study to determine how the behavioral characteristics of 
these two parameters resemble each other in synoptic time 
scales . 
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Several spectral studies of tropical cloud brightness 
[Tanaka and Ryuguji (1971, 1973), Murakami and Ho (1972a, 
1972b), Wallace (1971), Sikdar et al (1972), and Wallace and 
Chang (1972)] have been carried out previously. While some 
of these spectral studies yielded interesting results, 

Wallace and L. Chang (1972) suggested that the number of 
spectral peaks revealed by these studies are so numerous 
that, invariably, it would be extremely difficult to extract 
significantly useful information from such results. They 
also pointed out that the power spectra of brightness are 
highly variable in both space and time, and even neighboring 
grid points may exhibit quite different characteristics. 
However, the main purpose of the spectral analysis in this 
study is not to identify tropical wave disturbances as had 
been done in other studies. The purpose in this case is 
simply to use the spectral analysis technique as a convenient 
tool to decompose both brightness and divergence into their 
various time scales that are of synoptic interest in the 
period of this investigation. In addition, correlation 
analysis is applied, mainly to determine the daily spatial 
relationship, as well as to complement the cross-spectral 
analysis of the time series of both fields. 
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II. DATA AND METHODS OF ANALYSIS 



A. DATA 

The brightness and upper-level wind data used in this 
study were provided by the National Center for Atmospheric 
Research (NCAR) . The brightness data, acquired by NCAR from 
the National Environmental Satellite Service, are digital 
values representing variations from total absorption (0) to 
total reflection (10) of solar radiation in the visual 
channel. These digitized values are averaged for 5° x 5° 
latitude-longitude squares and adjusted for instrumental and 
calibration errors (Gray and Oort, 1974). The upper-level 
wind components (u, v) , which were used to generate the diver- 
gence time series, were acquired by NCAR from the National 
Meteorological Center tropical grid analyses. These winds 
are given at grid points approximately coinciding with those 
for the brightness data, and are the result of analyses based 
on all available rawinsonde records, aircraft observations, 
and a few data deduced from cloud drift as seen from satel- 
lite photographs. The 0000 GMT wind data are considered to 
be closest to (within 2 to 3 hours of) the satellite observa- 
tion time. 

The divergence field for the upper levels was derived 
directly from the horizontal wind components without any ad- 
justment by vertical boundary constraints. Maximum divergence 
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was found at the 200-mb level, consistent with the previous 
divergence prof iles found by Wallace (1971), Reed and Recker 
(1971), Williams and Gray (1973) and others. In these pro- 
files, the levels of maximum divergence are near 200 mb, 
indicating that this is generally the outflow level of tropi- 
cal motion systems associated with cloud cluster activities. 

Missing data were filled by a combination of simple and 
quadratic interpolation in time and a 4-point averaging in 
space. The projected values were restricted within the maxi- 
mum and the minimum of the observed field. 

The general area selected for this study is the tropical 
western North Pacific from the Equator to 25N and from 125E to 
175E. There are two reasons for this choice. Firstly, 
previous studies indicated that this is an area of strong 
cloud cluster activity. Secondly, the analyzed wind data 
in this region are more reliable than those in other tropical- 
locations because of the relatively denser rawinsonde network. 
Figure 1 shows this observational network in the region of 
investigation. 

The period of observation is approximately eight months, 

19 April - 20 December 1971. This period was chosen because 
daily records of IR radiation are available for the first 83 
days, namely: 19 April - 10 July 1971; hence, further investi- 
gation using I R- modified brightness could be pursued to shed 
additional light on the findings of this study. This is 
possible because, as mentioned in Section T, cloud-top 
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temperatures and cloud-top heights or the vertical extent 
of cumulus development can be estimated from the known IR 
field. Such estimates may then be utilized to modify the 
brightness pattern and, consequently, improve its correla- 
tion with that of the upper-level divergence or vertical 
motion. 

B. SPECTRAL ANALYSIS 

The BMD02T portion of the UCLA Biomedical Statistical 
Program Package was used for power spectral (and the asso- 
ciated cross-spectral) analysis. Both the brightness and 
divergence time series were first subjected to a high-pass 
filter which is described in Part D of this Section. A 20- 
day lag was used, which gives a frequency bandwidth of 
0.05 cycle day ^ (cpd) and a frequency limit of 0.5 cpd 
[Bendat and Piersol (1971), Panofsky and Brier (1965), 

Jones (1964), and Blackman and Tukey (1958)]. Spectral in- 
formation is then given at a frequency interval of 0.025 cpd 

Coherence squares, as listed in Tables I and II, repre- 
sent averaged values within each of the two dominant fre- 
quency bands, namely: 0.075 - 0.125 cpd (8.0 - 13.3 days) 
and 0.175 - 0.225 cpd (4.4 - 5.7 days), hereafter referred 
to simply as the 10-day and 5-day period bands, respectively 
The phase differences were realistically chosen to be repre- 
sentative (i.e., consistent phase differences produce a 
nearly average phase estimate, as used in previous studies) 
of each band in question. 
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To analyze the coherence square results, the question 
of statistical significance was resolved by taking into 
account that averaging the spectral estimates over three 
consecutive frequencies would amount to doubling the number 
of degrees of freedom implied by the individual spectral 
calculations. Consequently, in estimating the significant 
values at the levels indicated, the coherence square proba- 
bility distribution as compiled by Amos and Koopmans (1963) 
was used. 

C. CORRELATION ANALYSIS 

To complement the spectral analysis technique, correla- 
tion analysis was conducted, using the BMD02D of the afore- 
mentioned Statistical Program Package. The unfiltered data 
(as space and time series) of brightness and 200-mb divergence 
were subjected to correlation analysis. 

The space series cross-correlation analysis involved 66 
pairs of observations each day; this was done on a daily 
basis from May to October, which are considered as the summer 
period in the tropics. The time series, on the other hand, 
consisted of 246 pairs of observations for each of the 66 
grid points. In both cases, the significant correlation co- 
efficients, as listed in Tables III and V, were determined 
following Panofsky and Brier (1965). 

D. SPECTRAL FILTER 

As cited earlier, before conducting the spectrum and 
cross-spectrum analyses, a high-pass filter is applied 
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following Holloway (1958). The filtering technique is quite 
similar to those utilized in other investigations [Murakami 
and Ho (1972), Yanai and Murakami (1970), Yanai et al (1968), 
and Maruyama (1968)]. A normalized Gaussian weighting func- 
tion W with a standard deviation 6 of 25/3 days is used, 
n 

namely : 

exp (-t^/26^)dt , (1) 

where dt = 1 day , 

t = ndt , and n=0 ...+24 . 

The frequency response of this weighting function effec- 
tively filters out (89% efficient) periods beyond 36 or 25 
days (Holloway, 1958). As seen in Figure 2, this spectral 
response is ideal for frequencies beyond 0.075 cpd or, equi- 
valently, for periods not exceeding 13.3 days. Mathematically, 
this frequency response R(f) is expressed as a function of 
the standard deviation 6 and frequency f : 

R(f) = 1 - R^(f) , (2) 

where 

R^(f) = exp (“27T^6^f ^) . (2a) 

A weighted average is then calculated for every 49 (ap- 
proximately 66 ) consecutive data of the original time 

series of 246 days, thus resulting in a smoothed time 
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filtered data is then generated by the deviation of 

the smoothed time series from the observed time series, viz.: 
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III. RESULTS AND DISCUSSION 



Before discussing the results of this study, it might be 
useful to realize that the techniques of spectral and corre- 
lation analyses can be interpreted as some means for under- 
standing better those meteorological phenomena the physical 
nature of which are partially understood. Moreover, the re- 
sults of these techniques should be considered not purely on 
a statistical basis but also on the premise that they are 
compatible with synoptic and dynamic considerations [Yanai 
and Murakami (1970) and Panofsky and Brier (1965)]. 

A. SPECTRAL AND CROSS-SPECTRAL ANALYSES 
1 . Power Spectra 

Figures 3 and 4 show the power spectra at each of 
the 66 grid points for brightness and 200 mb divergence 
(div^QQ) , respectively, as generated by (4). An examination 
of these spectra indicates that the two frequency ranges: 
0.075 - 0.125 cpd (10-day period band) and 0.175 - 0.225 cpd 
(5-day period band), contain substantial segments of the 
total variances of each of the two series at most points. 
This means that during the period of this study, the diver- 
gence pattern at 200 mb associated with synoptic wave dis- 
turbances has characteristic time scales of 5 and 10 days. 
Hence, it is convenient to choose these two bands for 
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comparing the behavior of brightness and ^^^200 * Although 

the main concern of this study is not to define the spectral 
peaks, it should be recognized, nonetheless, that most of 
the major peaks do fall into these two bands. In this con- 
text, it might be significant to recall that previous 
studies [Sikdar et al (1972), Murakami and Ho (1972), Tanaka 
and Ryuguji (1973) and Gruber (1974)] found similar 2-band 
structures in the spectra of tropical cloudiness. It is 
apparent then that while there is a high variability in the 
brightness spectra, simple interpretation of a 2-band distri- 
bution is quite consistent with several studies and should 
not be overlooked. Additionally, one can notice a tendency 
for a 3-day spectral peak to appear occasionally; however, 
the peaks seldom appear concurrently in both parameters and, 
further, the fraction of the total variance involved is 
relatively small in most cases. 

2 . Coherency and Phase 

The coherence squares and the corresponding phase 
differences between brightness and ^^^200 each grid point, 

within each dominant frequency band, are listed in Tables I 
and II. The averaging process results in a degree of freedom 
of approximately 25, implying a significant coherence square 
of 0.17 at the 99% confidence level, 0.12 at the 95% level, 
and 0.09 at the 90% level. Only those values which are 
significant at the 90% level or better are tabulated. For 
the 10-day band (Table I), 24 grid points have significant 
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coherence squares, most of which indicate phase differences 
from in-phase to nearly 1/4 cycle between the two parameters. 
For the 5-day band (Table II), 20 points have significant 
coherence square values and, except for three points [(5N, 
135E) , (20N, 145E), and (20N, 150E) ] , there seems to be a 
similar relationship as in the 10-day band, although the 
larger (close to 1/4 cycle) phase difference is more common 
in the 5-day band* A relative comparison of these results 
for both frequency bands with those of Wallace (1971) for 
the 4- to 5-day waves indicate that, although the two para- 
meters are much more in-phase than out-of-phase in this 
study, the. coherence squares are somewhat lower and the 
phase differences are slightly larger. The reason for this 
difference could be explained in the light of the following: 
(a) smoothed, analyzed winds were used for computing diver- 
gence in this study compared with the direct radiosonde 
data used by Wallace (1971), (b) Wallace used a vertical 

motion field which was kinematically calculated from a verti- 
cally adjusted divergence field, as compar ed - wi th the non- 
adjusted divergence field in this investigation, and (c) the 
time difference between the two fields was ignored in this 
s t udy . 

3 . Horizontal Structure 

Another indication of correspondence between the two 
parameters in the 10-day and 5-day bands may be found by com- 
paring the horizontal structures deduced from i n t c r- lo ng i t ud c 
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cross-spectra for both parameters at each latitude, as shown 
in Figures 5-6, For each latitude, the 150E series is used 
as the base series to cross with other series at the same lati- 
tude. Each phase difference is plotted against the longi- 
tudinal distance from the base series and, unless the points 
are widely scattered, a best fitting line is drawn for each 
latitude weighted by the confidence level of the coherence 
square associated with each phase value. This technique has 
been commonly used in spectral wave studies for determining 
wavelengths and propagating directions. 

For the 10-day band (Figure 5), a westward movement 
in the brightness data with an apparent "wavelength" of 
about 50° (5600 km) could be inferred at all latitudes except 

the Equator, suggesting a "phase speed" of nearly 5° day ^ . 

This simple interpretation, however, seems to be inapplicable 
in the eastern portion of ION and 25N, and the easternmost 
part of 20N. Further, one may take note that the brightness 
field is more organized than the divergence field where the 
brightness results are best, viz.: 5N, 15N, and 20N. It 

seems significant to note, moveover, that another similarity 
can be found at ION, where an in-phase relationship is ob- 
served for the divergence series between 140E and 175E, and 
the brightness series between 145E and 170E. These results 
may be indicative of an east-west orientation and suggest 
the presence of an In ter tropical Convergence Zone (ITCZ) with 
intensity fluctuating around a period of approximately 10 
days. 



24 



The structures of the 5-day band (Figure 6) seem to 
be less regular compared to those of the 10-day band, with 
the divergence, again, less organized than the brightness. 

One may take not also that a partial best fitting line, at 
least, could be drawn for brightness at all latitudes, 
although the ’’wavelength** ranges from about 25® (2800 km) 

at the northern latitudes to nearly 50® (5600 kra) at the 

southern latitudes, thus indicating a **phase speed” range of 
approximately 5® to 10® day ^ . It might be of interest to 

note further that only two organized patterns, at 5N and 
20N, may be observed for the divergence but, in each case, 
agree well with those of their respective brightness fields. 

The fact that the brightness pattern is more organized than 
the divergence field, as revealed in the inter-longitude 
cross-spectra indicated in Figures 5 and 6, reaffirms our 
belief that brightness is qualitatively better than divergence. 

A similar attempt was made to make reasonable infer- 
ences from the inter-latitude cross-spectra, however, no 
significant in terpeta tion could be deduced because of the 
low coherences found in most cases. 

B. CORRELATION ANALYSIS 

1 . Correlation of Space Series 

The daily correlation coefficients of brightness 
versus div^^Q as a function of space (unfiltered data) is 
listed in Table III for the tropical summer, i.e.. May - 
October 1971. The sample size is 66, implying a signficant 
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correlation estimate of ^ 0.24 for the 95% confidence level, 
and ^ 0.20 for the 90% level (Panofsky and Brier, 1965). It 

4 

is interesting to note that of the 167 days with useful data , 
56 days (33.5%) indicate positive correlation at the 90% 
level or better while nine days (5.4%) tend to show a nega- 
tive correspondence. These results which, by and large, are 
similar to the findings for July 1969 (Edwards, 1973) suggest 
a weak positive correlation between the two fields on a daily 
basis^. One may notice, likewise, that for the "wet’* months 
(July - October), of the 109 days with useful data, 42 days 
(38.5%) show positive correlation while 5 days (4.6%) indicate 
negative correlation. On the other hand, out of 58 useful 
data for the relatively "dry" months of May and June, 14 days 
(24%) show positive correlation and 4 days (6.9%) indicate 
negative correlation. These comparative results appear quite 
consistent with the prevailing theories on the correspondence 
between precipitation and areas of low-level convergence or 
upper-level divergence. [See for example: Holton (1971) and 
Palmen and Newton (1969).] It should be realized, however, 
that these correlation coefficients are somewhat low as com- 
pared to those (inferred from the coherence squares) observed 



4 

Correlation estimates for which brightness and/or wind 
components (used to derive the divergence field) were 
originally missing were discarded. 

^Because of the many uncertainties in the computed diver- 
gence, tlie positive correlation listed in Table III (and 
Table V discussed in the next Section) should not be viewed 
as equivalent to the probability of finding a positive corre- 
lation between brightness and the actual 200- mb divergence at 
a given time (or space). 
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by Wallace (1971) for the 4- to 5-day waves at the KEP triangle, 
although they do indicate a tendency of in-phase relationship 
between the two parameters on a daily basis. Besides the 
possible errors involved in the kinematic computation of 
divergence and the time difference between the two fields, it 
is evident that the correspondence between the two fields is 
far from perfect. 

Moreover, it may be expected that the days of positive 
spatial correlation is partially due to the presence of strong 
convective systems. Table IV is a comparison of the frequency 
of positive correlation (Rf) between brightness data and 
200-mb divergence with the frequency of tropical storms (Sf) 
in the region of study from May - October 1971, which shows 
a nearly constant Rf to Sf ratio. However, it must be realized 
that usually tropical storms occupy only one grid point 
(horizontal scale - 5°) at any one time. 

2 . Correlation of Time Series 

The correlation coefficients between the time series 
of the two parameters under study (Table V) were similarly 
computed at all the 66 grid points in order to examine their 
spatial distribution. The sample size is 246 (unfiltered 
data covering the period 19 April - 20 December 1971), result- 
ing in a significant correlation estimate of ^ 0.12 at the 
95% confidence level, and ^ 0.10 at the 90% level. For the 
entire region, 25 points (37.9%) indicate positive correla- 
tion at the 90% level or better, while 3 points (4.5%) show 
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an apparent negative correlation. These results appear to 
be compatible with the daily computations discussed earlier 
and, likewise, are comparatively lower than the results found 
by Wallace (1971). We do notice, however, that the three 
negative correlations are in the southeast corner of the 
region and that all the interior points indicate positive 
relationship. Edwards (1973) examined the relationship be- 
tween brightness and vertical motion for the July 1969 data 
and found that the correspondence between the two fields is 
quite irregular in the eastern portion of the region. He 
suggested that this may be mainly due to the relatively weak 
cloud cluster activities in that area. Thus, if the region 
of study were narrowed to the area bounded by 135E and 160E, 
from the Equator to 20N, the daily correlations would be 
much higher . 
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IV, SUMMARY AND CONCLUSION 



In summary, this study is an attempt to investigate the 
temporal and spatial relationships between digitized bright- 
ness data and 200-mb divergence, as derived from the analyzed 
wind field, over a broad tropical region, namely: from the 

Equator to 25N, and from 125E to 175E. Spectral analysis was 
utilized to determine the dominant synoptic time scales for 
both parameters during the period 19 April - 20 December 
1971. To complement this technique, correlation analysis was 
applied at all the 66 grid points between the time series of 
each parameter for this 8 month period, and on a daily basis 
during the tropical summer: May - October 1971. 

Two common period bands, one centered near 10 days and 
the other near 5 days, have been observed. The cross-spectra 
between the digitized brightness and 200-mb divergence indicate 
a generally in-phase (< 1/4 cycle) relationship in these two 
bands, and further that the coherence squares are somewhat 
lower and the corresponding phase differences are slightly 
higher than those found by Wallace (1971) for the 4- to 5-day 
waves using direct radiosonde data. The horizontal structures 
of both parameters, as revealed by the in t e r- Ion gi t ud e cross 
spectra in the two bands, suggest that the brightness pattern 
is much more organized than that of the divergence. It is 
apparent, however, that the two fields resemble each other 



29 



whenever an organized pattern in the divergence can be found. 
Moreover, the zonal structures are less regular in the 5-day 
band than in the 10-day band. In addition, there is an indi- 
cation of a **waveleng th" of approximately 5600 km for the 10- 
day band and 2800 to 5600 km for the 5-day band, or a **phase 
speed** of about 5° day ^ and 5° to 10° day respectively. 

Correlation results for the space series suggest a relatively 
weak positive relationship between the two fields on a daily 
basis, and are consistent with the results for the time 
series. A somewhat higher correlation is found during the 
wet months of July - October than in the relatively dry months 
of May and June. 

The correspondence between the current brightness and 
divergence data as indicated in both the spectral and corre- 
lation analyses is far from perfect. However, in their 
entirety, the various results of this study do form a picture 
of positive relationship between the two fields. Thus the 
prospects of estimating the large-scale convection (vertical 
motion or divergence) field over a broad tropical region by 
the use of satellite data, as proposed by Wallace (1971), 
appear bright. Such estimates are believed to work best in 
areas of pronounced convective activity because of the 
expected influence of the presence of strong convective sys- 
tems. Moreover, until the present brightness data can be 
improved upon they are more useful in tropical wave studies 
than direct operational applications. The statistical results 
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of this study tend to support this point of view. For ex- 
ample, it was shown that the digitized brightness data 
reflect quite well the upper-level divergence pattern asso- 
ciated with synoptic-scale disturbances. Hence, it may be 
used to investigate the upper-level heating patterns since, 
by thermodynamic considerations, divergence is mainly 
related to heating in the tropical upper troposphere. 

Evidently the divergence data (as calculated from the 
smoothed, analyzed winds) contained several errors and, thus, 
might have reduced the correlation between brightness and 
divergence in this study. However, obtaining better diver- 
gence values from wind reports presently available in the 
tropics is conceivably difficult unless the observational 
network in the area is improved. This is precisely the very 
reason for current attempts to exploit the maximum utilization 
of satellite data. Also, problems do arise in utilizing 
brightness data to represent heating due to cloud cluster 
convections. Hopefully, these problems may be lessened by 
the availability of IR data which may be used to enhance 
the brightness data. The reason is primarily that the ver- 
tical extent of cumulus development can be determined from 
the IR information. Thus, an investigation on the relation- 
ship between IR-modified brightness and divergence is 
r e commend ed . 
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Same as Table I except for the 5-day band. 
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Day 

1 

2 

3 

A 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 



TABLE III 



Daily correlation coefficients (in hundredths) of satellite 
cloud brightness vs 200-mb divergence for May-October 1971. 
Values in parentheses are significant between the 90% and 
95% confidence levels, and values without parentheses, at 
the 95% level or better. The letter "M" represents missing 
data for either or both the brightness and the wind field. 



May 


Jun 


Jul 


Ana 




Oct 


27 


— 


(-21) 


M 


— 


- 


— 


24 


- 


28 


- 


- 


- 


- 


(-22) 


(22) 


- 


- 


26 


- 


- 


29 


- 


26 




- 


- 


28 


26 


H 


- 


- 


- 


- 


- 


- 


35 


- 


(22) 


- 


M 


- 


- 


-34 


28 


M 


42 


33 


- 


- 


(22) 


-31 


(20) 


- 


(21) 


- 


- 


- 


34 


36 


- 


38 


32 


- 


- 


34 


(22) 


45 


38 


- 


- 


(20) 


- 


32 


30 


- 


M 


M 


- 


- 


34 


- 


- 


M 


33 


- 


- 


- 


44 


M 


- 


- 


- 


- 


52 


M 


- 


- 


- 


- 


42 


- 


- 


-31 


- 


- 


(23) 


- 


- 


- 


— 


- 


- 


— 


— 


- 


- 


- 


36 


- 


- 


(22) 


(23) 


(20) 


(25) 


- 


- 


- 


40 


29 


- 


- 


(20) 


- 


41 


- 


- 


24 


- 


M 


(21) 


- 


- 


- 


M 


- 


- 


- 


(22) 


M 


- 


- 


-29 


(23) 


(21) 


36 


36 


-33 


- 


- 


- 


(20) 


- 


(-21) 


- 


- 


(20) 


- 


- 


- 


29 


- 


(-22) 


- 


M 


24 


39 

37 


— 


— 


— 
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TABLE IV 



Comparison of the frequency of positive corre- 
lation (Rf) between brightness data and 200-mb 
divergence with the frequency of tropical 
storms*(Sf) in the region of study from May to 
October 1971. 



Rf 

Sf 

Rf/Sf 



May 

8 

4 

2.0 



Jun 

6 

2 

3.0 



Jul 

14 

8 

1.75 



Aug 

7 

4 

1.75 



Sep 

13 

6 

2.2 



Oct 

8 

4 

2.0 



Average Rf/Sf = 2.1 



1. The term "tropical storms" in this table includes typhoons. 

2. Source: U.S. Fleet Weather Central/Joint Typhoon Warning 

Center, 1972: Annual Typhoon Report 1972 . Guam, Mariana Islands. 
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Correlation coefficients (in hundredths) of cloud brightness series 
vs 200-mb divergence series for 19 April - 20 December 1971. Values 
shown in parentheses are significant between the 90% and 95% confi- 
dence levels; values shown without parentheses are significant at 
the 95% level or better. 
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Figure 2, Spectral response function for the high-pass 
filter used in this study. 
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Figure 3. Power spectra of digitized satellite brightness data at the 
5° X 5*^ grid points. The ordinate is variance per unit frequency inter- 
val (B^ per 2 tt X 40 ^ day“^ x 10“^, where B is brightness unit with a 
range from 0 to 10) and the abscissa is period (days) . 
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Figure 3. (continued) 
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Figure 4. Same as Figure 3 except for 200-mb divergence and the unit 
for the ordinate is sec~* per 2ir x 40~*day~^x 10 
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Figure 5, Inter longitude 
cross-spectra of (a) bright- 
ness and (b) divergence for 
the 10-day band. The 
abscissa is longitude of the 
series that is crossed with 
the base series (150E) ; the 
ordinate is phase difference 
with positive values indicat- 
ing that the base series 
leads the other series. 
Plotting symbols indicate 
the significant coherence 
square at a given confidence 
level associated with each 
phase difference in the fol- 
lowing way: Blackened circle, 
^ 99%; circle with a cross, 
95-98%; open circle, 90-94%; 
small dot, < 90%. 



(a) 



(b) 
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Figure 6. Same as Figure 5 except for the 5-day band 
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